NEW  YORK  L'NtV'FP.StTY 
AFCRC-TN-57-760  WSTTimE  CF  f.  :  ICM 

ASTIA  DOCUMENT  No.  AD   133670 

la  W«M»V  Piac^  New  Yofi  3,  N,  Y. 

^>,     ^^    ^^      NEW   YORK   UNIVERSITY 


S    I    |J   |Y     73      Institute  of  Mathematical  Sciences 
:^      ^*-^      j<      Division  of  Electromognetic  Research 


RESEARCH    REPORT    No.    EM-107 


The  Inverse  Laplace  Transform 
of  an  Exponential  Function 


F.  M.  RAGAB 


CONTRACT  NO.  AFl 9(604)1717 
JULY,    1957 


NEW  YORK  UNIVERSITY 
Institute  of  Mathematical  Sciences 
Division  of  Electromagnetic  Research 

Research  Report  No,  EM- 107 


THE  INVERSE  LAPLACE  TRANSFORM  OF  AN  EXPONENTIAL  FUNCTION 

/ 
/ 

F.  rf.  Ragab 


f.Vi,   Ragab 

Morris  Kline 
Project  Director 

July,  1957 


The  research  repojrted  in  this  article  has  been  sponsored  by 
the  Air  Force  Cambridge  Research  Center,  Air  Research  and  De- 
velopment Command,  under  Contract  No.  AF19(60li)1717» 


Table  of  Contents 

Page 

1«   Introduction  1 

2«   Derivation  of  the  original  function  2 

3«   Derivation  of  a  subsidiary  formula  8 

h»       Darivation  of  the  case  m  «  2  9 

5b   The  asymptotic  behavior  of  the  original  function  10 
for  t  near  0 

Appendix  15 

References  16 
Distribution  List 


Abstract 


1  1 


n  mm 

The  inverse  Laplace  transform  of  the  image  function  p   Gxp(-a  p  ), 
which  arises  in  electromagnetic  diffraction  theory,  is  obtained  and  the 
asymptotic  behavior  of  the  original  function  for  values  of  the  independent 
variable  near  0  is  then  determined. 
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1,   Introduction 

In  electromagnetic  problems  one  often  has  occasion  to  study  pulse  fields, 
that  is,  fields  created  by  sources  which  have  a  Heaviside  unit  function  behavior 
or  a  delta- function  behavior  in  the  time  variable  as  well  as  fields  created  by 
the  same  sources  and  having  a  harmonic  time  behavior.  The  latter  fields  also 
prove  to  be  harmonic  in  time  (apart  from  a  transient)  and  these  time  harmonic 
fields  are  the  Laplace  transforms  of  the  corresponding  pulse  fields'-  ',  It  is  also 
true  that  if  one  expands  the  pulse  solution  of  Maxwell's  equations  (or  the  scalar 
wave  equation)  in  powers  of  (t-t  )  where  t  is  a  singularity  of  the  pulse  solution 
and  if  one  obtains  an  asymptotic  series  representation  of  the  spatial  behavior  of 
the  corresponding  time  harmonic  field  in  powers  of  1/co  where  co  is  the  circular 
frequency  of  the  source,  then  the  individual  terms  of  the  asymptotic  series  are 
essentially  Laplace  transforms  of  the  corresponding  terms  in  the  series  of  powers 
of  t-t^  (see  eq.  (U.S)  of  [8]). 

To  extend  the  knowledge  of  the  relationship  between  time  harmonic  fields  and 
pulse  fields  and  between  the  asymptotic  series  described  above  and  the  correspond- 
ing series  of  powers  to  fields  diffracted  by  smooth  bodies,  it  now  appears 
desirable  to  utilize  some  Laplace  transforms  which  have  not  been  established. 
The  present  paper  is  concerned  with  determining  the  original  function  whose 
Laplace  transform  is 


(1)  p    expl-a  '  p  ^   , 


where  p,  as  usual,  is  complex,  R(p)  >  0,  R(a)  >  0,  R(a)  >  0  and  m  »  2,3,h,...  • 
It  is  expected  that  the  result  will  be  useful  in  electromagnetic  diffraction 
problems  but  the  present  paper  will  be  devoted  to  the  purely  mathematical  problem 
of  determining  the  original  fvinction  and  also  the  asymptotic  behavior  of  the 
original  function  for  t  near  0» 

When  m  =  2  the  function  (1)  is,  of  course. 
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(2)  p-lexpp/2pl/2j 


and  the  original  function  is  known  (see    [l] ,  p.  2U6,  eq.  9).  It  is  in  fact  the 
function 

i-«-^  -*  111 

(3)  2^      '^    „    ^t-e     ^    ^2a-lC2"  ^  -^   t"  0 


where  R(a)  >  0,  R(a)  >  0  and  T)^{z)  is  the  parabolic  cylinder  function  (see   [?] ,  p.   392) 

^+  ^       ^  2 

D^(z)-2^*^z"^   W^       1         l(|-) 

(U)  o-v-l  V  2/,  -  2 

The  function  appearing  on  the  right  side  of  this  last  equation  is  MacRobert's 
E-function  whose  definition  and  properties  are  to  be  found  in  the  literature 
(see  [3],  p.  3U8  and  352  or  [6],  pp.  203-206)  and  about  iriiich  more  will  be  said 
below.  Hence  the  major  result  of  this  paper  is  a  generalization  of  the  Laplace 
transform  relationship  between  the  functions  (3)  and  (2)  to  the  case  where  the 
exponents  in  the  exponential  part  of  (2)  are  l/m.  The  result  will  be  stated  and 
derived  in  Section  2.  After  deriving  a  subsidiary  formula  in  Section  2  we  show 
in  Section  k   that  the  general  result  of  Section  2  reduces  to  the  special  case 
when  ra  ■  2.  In  Section  5  we  derive  the  asymptotic  behavior  as  t  approaches  0  of 
the  inverse  Laplace  transform  of  the  function  (1). 

2*     Derivation  of  the  original  function 

We  present  first  for  convenient  reference  several  facts  about  the  E-function 
and  other  relationships  which  will  be  utilized  later.  We  begin  with  some  facts 
about  the  generalized  hypergeometric  function. 

The  symbols  (k|0)  and  k|n),  where  k  may  be  real  or  complex,  are  defined  as 
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(kjO)  -  1 

(kjn)  -  k(k+l)(k+2)...(k+n-l),       n  -  1,2,3, 


•  •  •  • 


Than  the  generalized  hypergeometric  function  FCpja  ;q}p  :z)  is  defined  aa 

C  a 

00   (cL,}n)(a5;n),..(a  ;n) 
F(p}a  :qiP  :z)  -  T"   —i, ij r—^ r-    z" 


wherein  it  is  assumed  that  none  of  the  p's  is  zero  or  a  negative  integer.  An 
alternative  notation  for  this  F-function  is 

When  p  »  1  and  q  ■  1,  the  generalized  typergeometric  function  F(l}a:l}o:z)  is 
often  denoted  as  F(aiP}z)  or  as  ,F,(a}p;z), 

If  p  <  q  then  the  E-function  is  defined  as 

r(a^)r(a2)...r(a_)  1 

E(pia  :q;p  :z)  -  — — ^  F(p}a  :q;p  :  -  -  ) 

^    «      r(Pi)r(P2)...r(Pq)      ^         ' 

where  z  /  0,  If  p  >  q+l-  and  this  is  the  casa  which  will  be  of  interest  here, 
and  When  |arap  z|  <  n,  then  the  E-function  (see  p.  353)  can  be  shown  to  be  K  * 

E(pja  :q;p  :z)  -  f   fr'  r(as-a^)<  f   P^Pt'^r^  ^(0^)  ♦ 
r«l  s"l  t=l 

z  ^F  1 

°r'°l*-^»***  *  •••»  °r'"^p*^ 

The  definition  of  the  F-function  for  this  case  is  given  in  the  Appendix  where 
it  is  utilized. 


-  u  - 

where  the  dash  in  the  product  sign  signifies  that  the  factor  for  which  s  ■»  r 
is  omitted  and  the  asterisk  in  the  F  function  means  that  the  parameter  a  -a  +1 
is  omitted,   so  that  there  are  only  p-1  parameters  in  the  lower  set.     When  there 
are  no  p  's  in  the  E-function,  that  is,  when  q  «  0,  then  equation  (5)   holds  under 

3 

the  conditions  that  the  factor  in  the  braces  is  omitted  and  the  parameters 

a  -p.,+1,   ...,  a  -p  +1  are  omitted  in  the  F-function. 

To  familiarize  ourselves  with  the  E-function  the  follo^ying  relations  may  be 

worth  notinp  though  they  will  not  be  used.     From  the  definition  of  the  E-function 

for  the  case  p  <  q  and  from  equation  (5)  it  is  clear  that  the  E-function  is 

immediately  related  to  the  generalized  hypergeometric  function  and  reduces  to 

simple  expressions  in  the  ordinary  or  Gauss  h^ypergeometric  function  when  p  «  2 

and  q  «  1.     For  p  •  1  and  q  ■  1  it  is  ailso  evident  from  the  definition  that  the 

F-function  reduces  to  the  confluent  hypergeometric  function  or  Rummer's  function. 

The   case  p  -  2,  q  ■  0  yields  the  following  relation  (see    [3] »  p.  351,  eq.   (lU)) 

1  1  2 

cos  nnE(x'  +n,  «•  -n::2z)     -    /2nz     e     K  (z) 


where  K^(z)  is  the  modified  Bessel  function  of  the  second  kind.  It  is  also  true 
E(^  -k+m,  ^  -k-m::z)  -  f  (7  -k+m)r(2-  -k-m)z"'^e^  W^^  ^(z) 


n 
that 


^ere  W.      (z)  is  the  Whittaker  function  (see    (3],  p.  3$1,  «q.   (15)).     Finally  it 
is  immediate  from  the  definition  of  the  E-function  that  for  p  »  0  and  q  -  0  the 
F-function  is  just  e~       .     More  complicated  parameters  in  the  E-function  lead 
to  the  equivalence  of  the  E-function  with  products  of  Hankel  functions,  with 
Lommel  functions,  and  with  products  of  Whittaker  functions   (see    [3],  p.  395). 
We   shall  utilize  the  relationship 


(6) 


CO  , 

/        %     a     , -1 


^-X  ^  p+1     E(p;a^:q;Pg!z/X)dX     -     E(p+l;aj.:q}Pg:z) 


o 


where  R(a     ,)  >  0.     This  equation  will  be  proven  in  Appendix  I, 
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other  equations  we  shall  use  are 


mil 

-raz 


smnz 


(7)  r(z)r(2+  b  ...  r(z*  ^)  -  (2n)^  ^  m^    r(inz) 

A  In 

which  is  the  multiplication  equation  of  Gauss  and  Legendre  (see  [6],  p.  U,  eq,  (11)); 

(8)  r(z)r(i-z)  • 

(see    [6],  p.   3,  eq.   (6))   ; 

(9)  ^Fj(ajpjz)     •     e''  ^F^(p-a}p}-z) 
(see    [3]>  p.   3U6,  eq.   (5));and  the  conmon  trigonometirLc  identity 

(10)  sin  -  sin  —     .*.     sin  n     ■     2       m* 

mm  m 

We  shall  be  interested  in  the  E-function  when  there  are  no  p  's,  that  is, 

s 

1  2 

when  q  =  0  and  when  the  p  parameters  a     are  the  m  parairaeters  a,  a  ♦  — ,  a  +  —  , 

X^  IQ  lu 

•••.  a  +  •  Then  we  shall  use  the  accepted  notation 

m 

^/  1  2  m-l         X 

E(a,a+  -,  a*  -,    ...,  a+  -^  ::z)   , 

From  (5)  with  the  conventions  agreed  upon  in  connection  with  ($)  we  may  say  that 

E(a,a.  i,...,a*  ^  ::z)  -  £      TT  TC"  s^'TT       r(|)r(a*  |) 

r"0     S"l  t"l 

(11) 

r  r 

a*  —        /a  +  -  ;  z 

^       "  ^  (  !♦  E."  1*  III  !♦  1  1-  1  1-  i  i-»lL-^ 

wherein  we  have  split  the  first  product  on  the  right  side  of  (5)  into  two 
separate  products*  The  F-f unction  on  the  ripht  side  of  (11)  is  sometimes 
also  accompanied  by  subscripts  thus:  iF__t  * 

We  shall  now  establish  the  following  Laplace  transform  theorem  in  which 
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the  original  function  is  the  function  of  t  and  p  is  the  Laplace  transform 

parameter  whose  real  part  is  greater  than  0: 

■  e-f"  r   ^  Ua,a*   i  .  a*  | a*  i!^  ::  ^  )dt. 

i,-i  m  t 

(12)  "  1  ^  m    1  ,^„  ^       , 

,^  ,7  *  ?  -  7  -"°  a  a-1    f  l/m  1/m  1 
(2n)     n      a  p    exp  -a    p 

where  R(a)  >  0,  R(a)  >  0,  m  -  2,3,U,...,  and  Y       means  that  in  the  expression 


E. 


following  it  i  is  to  be  replaced  by  -i  and  the  two  expressions  are  to  be  added. 
To  prove  equation  Cl^)  we  shall  transi'orm  the  left  side  into  the  right.     We 

consider  therefore  the  left  side.     We  substitute  t  -  X/p  and* apply  (6)   in  which 
we  let  the  value  of  a     .   -  1.     This  yields 

,00  in 

{     -pt  ^ —  1  -/  12  m-1  ae       \.. 

e  *^     >  T-  E(a,a+  -,a+  -,..., a+  — --  ::  — — -  )dt     - 

J  .•—.  i  '       m'       m  m  m. 

We  now  express  each  E  function  on  the  right  side  in  terms  of  the  value  given  by 
(11)  and  combine  the  two  resulting  expressions  by  factoring  out  common  terms. 
The  terms  which  result  from  2*^*'"/"  in  (11)  are  transformed  by  means  of  (8).  We 
also  use  the  definition  of  the  generalized  hypergeometric  function  to  cancel 
factors  in  ntunerator  and  denominator.  Hence  the  right  side  of  the  preceding 
equation  becomes 

m-1 


^      m        r«o 


(2E)"  r(-  |)r(-  '-^)  -n-  i)r(i)r(|)-r(^) 

m 


(»a)  •  A.i(ii*  „^,i*  ^,....x.  i,  1-  i.  1-  |...A-  =is  ;^ 


) 

-1 
m 


where  the  notation  in  F  -  means 

0  m-1  ,. 


(-1)  ap_ 
T      1     2 


m 
m 


m'     m 


m-l-r 
m 
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Now  apart  from  the  power  (ap)  before  the  summation  in  (13a)  the  expression 

contains  powers  of  (ap)  where  s  =  0,1,2,,,,  in  the  generalized  hypergecmetric 

m-1  / 

function  each  multiplied  by  ^     (ap)  '    .     Let  us  however  concentrate  for  the 

r»o 
moment  on  a  particular  value  of  s  and  on  a  particular  value  of  r,  the  former  in 

the  range  0,1,2,,,,  and  the  latter  in  the  range  from  0  to  ra-1.  By  using  the 

series  meaning  of  F  ,  we  find  that  the  coefficient  of  (ap)  '""  *  in  (13a) 
o  in— X 

equals 

2n  ^-ma-r-mS(^p)a(_^)ras  . 

^^^^       r/p(-i)r(i+  i)--T(--  -)r(i*  -)|  {r(-)r(i-  i)-"r(^^2^^^)r(i-  ^!^^^^)\ 

V   m  '    m^   '  ^  m  '  ^   m  J  L  '  ^m  '  ^   m    '  ^  w   '  ^ m   J 

.'.  (r(i-  ^  .s)--r(i*  i'-s)]  (pd-  i^s)...r(i-  ^  *s)j 


1 


We  now  use  (8)  on  the  numerator  of  this  expression  to  obtain 

2n  -ma-r-ms/  \a,   - ^r+ms  n 

F  °       ^^P^  ^'^^     f  ,  Ti  .   ?n ,  rnl  r  ■  n  ■   2n .  m-r-1   1 

'^  \  sin  —  sin  —  •••sin  —  Usin  —  sm  —  •••sin n  V 

Lm  m  mj\ra  m  mJ 

(15) 

1 

P(i.3)[r(i-  —^  *s).,.r(i- 1  *s)}  [ra*  I  *s).„r(i*  i  -s)] 

We  now  use  the  simple  identity  sin  —-  »  sin(jt n)  -  sin n,  for  r  ■  1,2,,.., 

in         m         m  »»» 

m-1  in  the  first  denominator  and  rewrite  the  factors  in  the  second  denominator 

with  a  trivial  change  in  form.     Hence   (l5)  becomes 

m-1 
2tt    -ma-r-ms/     xO/  ,xr+ms n 

p-  »  (^P^   ^"^^  ■      n     ■      2n   ,.,    .     rn     ,      r>l — TT^ .      m-1 

'^  sin  —  sin  —  •••sin  —  sin  n  sin  n,,.8in  n 

mm  mm  m  n 

(16) 

1 

P(i* i+s)r(-  +  -+s).^^r(i  +  ^*  s *^ 

'    ^m       m  '    ^m       ro       '        '    ^m       m  m   ' 

Use  of  relations  (10)  and  (7)   converts   (16)  to 
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-  m  (ap)   (-1)  -^^     -T-- jp-T^ 

^  in2  «  -1-r-ms         7 '  5- 

J.  m^  (2n)^      '^    Pd+r+ms) 

—  +s 
Thus  the   coefficient  of   ''ap)'"         in  (13a)  becomes 

1      "I       1  ~^-,- 

?+  ■s-    ~?r  -ma  ^    /  1  \r*ins 

2       ?  a  G-1  (-1) 

Ir+ms 


will  oh  equals 

1  1  f  1  i    i 

i  ♦  y    -i-ma  ,  -  (r+ins)  _  m     m 

(2n)  m  a°p°"Mcoef,  ofCap)""  in  the  expansion  of  e         ^ 


If  we  now  consider  the  successive  powers  of  (ap)  '     in  the  series  expansion 
^l/m  l/n 
of  e'         P        we  see  that  we  shall  get  first  of  all  the  set  of  powers  from 

n  =  0  to  n  -  m-l  which  frives  the  powers  of  (ap)  for  s  -  0  and  r  from 

0  to  m-l.     The  powers  from  n  =  m  to  n  »  m  + in  the  exponential  series  gives 

T*  m 

the  powers  of  (ap)     for  s  «  1  and  r  from  0  to  m-l.  More  generally  every 

—  +s 
power  of  (ap)"    which  is  called  for  in  (13a)  (apart  from  the  factor  (ap)*^  in 

1/m  1/m 
front)  is  contained  in  the  series  for  e~^   ^    and  only  those  powers  are  con- 
tained in  the  exponential  series.  Thus  (13)  becomes 

e    >    T  F(a,a+  r  »a+ 


Ptr-  \.<.  ..\    ..l  ^^„^^.^  ,,^)  ^x,    - 


r   1  F(a,a.  i  ,c 


m  »'•-»-   m     m. 

n  \j 


1  +  m      1   n,« 

/ox  ,^  X?   ?  "  ?  ''"^  a  a-1    ,   1/m  l/m^ 

(16)  (2n)    ^  m      a  p    exp(-a  '  p  '  ), 


which  is  the  desired  result  as  stated  in  (12), 


3,   Derivation  of  a  subsidiary  equation 

We  wish  to  show  that  the  general  transform  (12)  or  (iB)  reduces  in  the  case 
m  =  2  to  the  Laplace  transform  relationship  between  (3)  and  (2).  For  this  purpose 
we  shall  first  prove  a  subsidiary  equation,  namely. 


'  9  ' 


5    .   a+p-1  -z  . 

(19)     £1115 2 E(l-a,l-p::z)  -  E(a,p:  ize^")   -  E(a,p:  :ze'^"), 

r(i-a)r(i-p) 

where  the  notation  in  the  F-functions  is  explained  by  (11). 

We  re-express  each  of  the  two  E-functions  on  the  ripht  side  of 
(19)  by  means  of  equation  (11),  add  similar  terms,  and  obtain 


E(a,p::ze^") 


E(a,p::ze'^") 


(20) 


2ni 


/Pi   -z 


rO-g)   ,a     p  f'^J  -^)    ^    Via-?)     F    f^'*   '"•  \ 


T 


We  now  use  relation   (9)  to  write 


2Ttiz 


a-M  .    f   r(a-p)r(l-a)zl-%F,(J3:') 


r(i-a)r(i-p) 


l-pjz 


by  applying  (11)  we  obtain  (19)  above.  Hence  the  subsidiary  equation 
is  established. 

U.   Derivation  of  the  case  in  °  2 

For  the  case  ui  =  2  the  general  result  (12)  tells  us  that  the  transform  of 


(21) 


J  E(a,a+  J  ''   -Tjt~  '  '   I  ^(g>2-^  J 


ae 


-in 


TTT 


is 


(22) 


3   1 
(2n)^2  "^ 


•2a 


1  1 

a  a-1    /   2  ?\ 
a  p    exp(-a  p  )  . 


We  shall  show  that  the  transform  relationship  between  expressions  (21)  and  (22) 
agrees  with  the  known  relationship  stated  in  expressions  (3)  and  (2). 
We  use  the  equality  (19)  to  transform  (21)  and  obtain 


wherein  we  let  x  =  «■  -  a  we  may  simplify  the  preceding  expression  to 
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2a-  »•  »  -2a 

We  now  use  equation  (U)  wherein  we  let  v  =  2a-l  to  obtain  from  the  preceding 
expression 

;„(g)'°'  ?  t'^  ''°     3/^,  „V^r(l-2,)   a/at.     /  a^/^   \ 
Dy  means  of  the  duplication  formula  for  the  gamma  function,  namely, 

r(2x)  -  r(x)r(x*  \)  - — 

X-  -  a  we  may  simpli; 

If  we  now  eliminate  constant  factors  in  this  expression  and  in  (2U)  we  see 
that  the  Laplace  transfonn  relationship  between  (21)  and  (22)  is  identical 
with  the  Laplace  transform  relationship  between  (3)  and  (2), 

5»   The  asymptotic  behavior  of  the  original  function  for  t  near  0 

It  is  known  for  a  wide  class  of  original  functions  and  their  Laplace  trans- 
forms that  the  asymptotic  behavior  of  the  original  function  f or  t  near  0  determines 
the  asymptotic  behavior  of  the  image  function  as  |p |  becomes  infinite  (see  [5], 
p. 137),  With  a  possible  extension  of  this  existing  theorem  in  mind  we  determine 
here  the  asymptotic  behavior  of  the  original  function 

-  i  -  "  J  +ma  -a  T-   1  i,/-     1     2        m-1    ae^"  v 
(25)      (2n)  ^   ^m^     ^   JI,^  T  «(«>-*  m>  °*  P  -  •>-*  —  ''  ^  >' 

We  utilize  a  known  integral  for  the  E-f unction  called  the  Barnes  integral 
(see  Ijl ,  p.  37U).  For  our  case  where  the  E-function  contains  no  p  's  and  the 
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.  12  m-l 

p  parameters  a     are  a,  a+-,a+-,    ....   a+  ,  then 

r  a      n  m 


-,/     1    2        m-l    ae^"v    1 
E(a,a+  -,  a*  -,  ...,  a+  -;:-  ::  —;^)   -  jj^ 


^  »-l     ,.     „  in  s 

V   \/ae 


n'    n'    '     n  m.         

m  t       / 


V"0  ni  t 


where  the  path  of  integration  may  be  taken  along  the  imaginary  axis  with  loops, 
if  necessary,  to  have  the  poles  a   +  v/m  lie  to  the  right  of  the  contour.  If  we 

consider  the  summation  ^   t  E  called  for  in  (25)  we  note  that  the  factor  e  "^ 

i,-i 
will  appear  in  the  integral  representation  of  the  first  term  of  the  sum  and 

"ilTS 

the  factor  e  will  appear  in  the  integral  representation  of  the  second  terra. 

Conversion  of  these  factors  to   oosine  and  sine  temis  permits  combination  of  the 
two  integrals  and  the  use  of  (8)   gives  the  result  that  (25)  equals  AI  where  I  is 
given 


1    r  i?ii  r(a+  ^-s) 


(27)        I  .  -  ^      TT       ^.V      (4-)'  ds 


THT 


v-o 


Pd-s)  m^t 


with  the  same  path  of  integration  as  before  and  where 

1ml 
A     »     (2n)  m  a       , 


We  show  next  that 


11  ^_  m 


(28)       exp|  (n,-l)(-i-)"^  I        (2n)  "  ^(m-D^^^  i  .  -  ^         S(8)(^)''d8 
m\  J  m  t 


where 

00        II       I    ^       m-l       m 


s(=)  ■  E    ^ 


u»o 


Tr©r(i.^-s) 


r-1 

when  R(s)  >  R(a-  ■*-) .     To  see  this  we  first  substitute  the  value  of  S(8)  in  (28) 
and   the  rif?ht  hand  side  becomes 
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(29) 


1 
7^ 


00 

L 

u»o 


m-1 


v»o 


m-r 


(4-)  ds   . 


r«l 


in  t 


We  interchange  integration  and  summation.     Then  (29)  becomes 


v«o 


m-l       m 


a 


K  r(^)r(i-^-) 


(-2-)     ds     . 

m. 
m  t 


r-l 


In  (30)  we  chanee  s  to  s  +  — -•  and  obtain 


00 


u 
a   \m-l  1 


m-l 

/  ¥  r(a*j-8) 


u-o  rat  -|y  p(liirj     y 


v«o 


m 


r-l 


r(i-s) 


(-J-)  ds. 

m. 
m  t 


ra- 


This  equals 


u 

a  \nP 


00 

u"0     m  X 


r(^'r(^)-r(^>S' 


I  . 


We  now  use  (?)  vjith  b-1  instead  of  n  to  obtain 


00 

E 

u«o 


m. 
m  t 


1   -,u 
5=T 


la 


(2n)    ^       (m-l) 


1 
w  +u 


r(u+i) 


I, 


which  equal '^ 


1- 


m 


00 


1     U 


(2.)    '^(m-l)^/^E    ^("-D^-T'"" 


u«o 


m  t 


¥ 


or 


1- 


m 


(2n)"  ^  Cffl-1)-'-/^  ejqp 


With  (28)  established,  we  utilise  a  known  asymptotic  expansion  of  the 
expression  on  the  right  side  of  (28),  which  holds  when  |   ^  |  is  large 
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(see   [Uj,    p.296-7)»     By  putting  all  the  factors  on  the  left  side,  of  (28),   except 
I  on  the  right,   we  conclude  that  for  |t|  near  0,   or  when  rrr  is  large, 


AI  /^    A  exp 


-(m-l)(-^) 

m. 

m  t 


1    -r 
a  \m-I 


1     (         '^\ 
r_a_-|iS=T^"'°"  ?^ 

'-ra  t-J 


(31) 


m 


(2n) 


-1 


\ 


M, 


farl 


a  \m~1 


m 
m  t 


rn  t 


+   ... 


wherein  the  coefficients  M_,  M^,    ...  do  not  depend  on  t.     Of  course  the  leading 
term  in  this  asymptotic  expansion,  in  view  of  the  value  of  A,  is 


_  1       »  +nia 
(2n)  -  '  a       exp 


/ra=T 


-(„,-l)(4.)^ 
in. 
n  t 


— m  t-l 


1     (  1^ 


We  now  consider  the  special  case  of  m  »  2  and  show  that  the  result  (31) 
agrees  for  this  case  with  a  known  result.  We  suppose  that  Ir-rl  is  large  which, 
of  course,  will  be  the  case  for  t  near  0,  We  substitute  m  »  2  in  (31)  and 
taking  account  of  the  value  of  A  obtain 


(32) 


2  a 


1  1 


+  2a 


-aAt(a.,"-?„„,?^?   'V".(^), 


where  J(pr)  tends  to  1  as  Ij-r-l  tends  to  infinity. 

This  result  agrees  with  the  known  asymptotic  expansion  of  the  original 
function  (3)  which  holds  for  the  case  n  «  2  ,  because  (see  [2],  p.  392,  line  2) 


.-'H^) 


°"7 


^2a-l(^/  ^  J  "  „V^p(i.2a)  "'^^^^  ^^^  "°*  ^"°  ''''  ^^    * 


1^^ 


However  (see  [3],  p.  3^1,  formula  (13)), 

E(i  -a,l-a  ::  A-)  -^  r(T-a)r(l-a)  J(i^) 


UT 


W^ 


-  lU  - 


where   JCt-t-)  — >  1  as  t  — >  0.     Hence 


{n^\  ^""'?t'" 


-2o,a  ,""? 


Use  of  the  Legendre  duplication  fonmila  for  the  gainma  function  gives 


Hence  the  asymptotic  value  of  the  original  function  (3)   is 

1  1 

o?      °    -1/2  .-a  ,a  x°  "  7      -a/Ut    ,/ a  x 

which  agrees  with  (12)« 
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Appendix 

We  now  prove  a  relationship  used  in  the  text  projjer,  namely, 
00  ^       _^ 

(1)  I  e"^  n  P*       E(pjo^:qjpg:3/ii)diJi  -  E(p+l;aj.:q}pg:2) 


when  R(a     ^)  >  0.  (This  relationship  is  stated  as  problem  106  on  p.  39li  of    f3]  .) 

Proof:     If  p  >  q+1,  we  use  the  definition  of  the  E-function  ([3],  p.  352, 
eq.   (20)),  namely, 


K(Pi.,=qiP,=.)  -  r(«p)  ]  J^r(p„-<.„) "  fr  j  x°"'  a-x/-'""'  dx„ 


o 

(2) 


o  00       .  .  ,00       ^  - 

p-2      /      -X       a  -1  /      -X     -     0.     ,-1  -a 

fr    (  «  "C  <^\^  '   S-1   (i*^L^...Vi/^^  Vi 

n-q+1  J^  i  V  y- 

where  R(a   )  >  0,  n  «■  l,2,,..,p-l,  R(p  -a  )  >  0,  n  ■  1,2, ...jq.     We   substitute 
this  expression  for  E  given  by  (2)  with  of  course  z  replaced  by  z/|x  into   (l). 
We  then  merely  place  the  terms  involving  |i  and  the  integration  with  respect  to 
\i  within  the  innermost  integral.     The  resulting  expression  is  seen  to  be,  by 
the  very  definition  of  the  E-function, 


(3) 


E(a2»°2»*'**°p-1'  %*V  °p"^^^6*^^   * 


However,  equation  (5)  of  the  text  shows  that  the  E-function  is  symmetric  in  the 

a     and  p  .     Hence  expression  (3)  is  equal  to  the  right  side  of  (1). 
r  s 

If  p  <  q,  we  use  as  the  value  of  the  E-function  the  definition  for  this 
case   ([3],  p.  352,  eq,   (21)),  namely. 


r(a^)r(a2)...r(a_)  1 

E(p;a   :qip   :z)     - ^       ■    F(pja   tqjp   :    -  -)  # 

'    ^        r(Pi)r(P2)--.r(Pq)      "^    ^    ' 

We  substitute  this  value  of  E  in  (1),  expand  the  generalized  hypergeometric 
function  into  the  series  which  defines  it  and  integrate  term  by  term.  In  view  of 
the  definition  of  the  Gamma  function,  the  result  follows  at  once. 
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